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This report is based on the work done together with J.B. Kogut and C.G. Strouthos. We study a Yukawa theory 
with spontaneous chiral symmetry breaking and with a large number N of fermions near the finite temperature 
phase transition. Critical properties in such a system can be described by the mean field theory very close to the 
transition point. We show that the width of the window of non-trivial scaling is suppressed by a certain power of 
1/N. Our Monte Carlo simulations confirm these analytical results. We discuss implications for the chiral phase 
transition in QCD. 



1. Introduction 

Global symmetries provide important informa- 
tion about the properties of quantum field theo- 
ries. In QCD with 2 massless quarks the global 
symmetries form the group: SUa(2) xSU(2)y 
xUb(1). This fact, together with the ideas of 
dimensional reduction and universality gives us a 
prediction for the long-scale behavior at the chiral 
symmetry restoration transition in QCD 

The idea of dimensional reduction is based on 
the observation that at the second order phase 
transition a certain correlation length in the sys- 
tem diverges and eventually becomes longer than 
the extent of the Euclidean time. The long- 
wavelength fluctuations are then squeezed to a 
d — 1 dimensional "pancake" . The universal- 
ity tells us that the long-wavelength behavior of 
a system depends on the global symmetry and 
the dimensionality of the system. In particular, 
for the QCD with 2 massless quarks one should 
find 0(4)=SU(2)xSU(2) exponents at the chiral 
phase transition [Q. 

Such an argument, based on universality, does 
not tell us how this universal critical behavior sets 
in. To answer this question one has to study the 
dynamics. In this work we study how this univer- 
sal behavior sets in for theories with large num- 
ber, N, of fermion species. 

The fermionic fields, unlike the bosonic, do 
not survive dimensional reduction. One way to 
see this is to recall that fermion fields have an- 
tiperiodic boundary conditions in the time di- 



rection, which suppress fluctuations when their 
correlation length exceeds 1/T. In another lan- 
guage, the effective d — 1-dimensional theory for 
the bosonic fields near T c is a classical statisti- 
cal theory. Fermion fields do not have a classical 
limit f§. 

However, the fermions may affect the way that 
the universal behavior sets in. Several authors 
observed (^^] that in the large N limit the ex- 
ponents of the finite T chiral transition in the 
Yukawa model are given by the mean field (MF) 
theory. Monte Carlo simulations confirmed this 
result 0. How does it reconcile with the univer- 
sality arguments? The answer is that the large 
N description has its region of applicability. The 
scaling for the correlation length in the window 
between 1/T and N x /T, x > 0, is indeed given 
by MF. However, when the correlation length ex- 
ceeds N x /T the universal scaling sets in. 

2. The model 

To derive the crossover exponent x we study || 
a Yukawa theory with large number N of fermion 
species. It can be thought of as an effective de- 
scription of degrees of freedom in QCD which par- 
ticipate in chiral symmetry breaking. We choose 
the symmetry group Z(2) instead of 0(4) to sim- 
plify the argument. The theory is defined by the 
Lagrangian: 
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and is regularized by some momentum cutoff, A. 
There are two other important scales in the the- 
ory: the temperature, T, and the mass m of the 
thermal excitations of the scalar field, m van- 
ishes at T c . Therefore, near the finite temperature 
phase transition we have the following hierarchy 
of scales: A > T > m. 

The renormalization group (RG) evolution of 
the quartic self-coupling of the scalar field, A, 
from the scale A to the scale T is governed by 
the RG equations of the d-dimensional quantum 
Yukawa model. After that, at the scale of T, we 
pass through a crossover region: the fermions and 
nonzero Matsubara modes of the scalar fields de- 
couple. The evolution below T is governed by 
the RG equations of the scalar (j> 4 theory in d — 1 
dimensions. 

If the window of scales between A and T is 
wide the renormalizcd coupling A at the scale T, 
A(T), is close to the infrared fixed point of the d- 
dimensional Yukawa theory. In the large- N limit: 

A(T) ~ ^ for 2<d<4. (2) 

The case d = 4 is special. The infrared fixed 
point is trivial and is approached logarithmically 
as A/T -► oo: A(T) ~ l/[Yln(A/T)]. For large 
N this coupling is small. As we shall see shortly, 
this is the reason why the non-MF critical region 
is reached only very close to the phase transition. 

The quantitative relation between the size of 
the non-MF critical region, the Ginzburg region, 
and certain parameters of a given system is known 
as the Ginzburg criterion. In superconductors 
such a parameter is a small ratio T/Ep, i.e., the 
width of the Ginzburg region is suppressed by a 
power of this parameter. In a field theory with 
large number of fermions, such as (0), such a pa- 
rameter is 1/N. 

The MF approximation breaks down because of 
self-inconsistency when the fluctuations become 
large. The size of the corrections to the MF is de- 
termined by the value of the effective self-coupling 
of the scalar field. Since the dimensionality of the 
coupling of the d — 1-dimensional scalar theory, 
Ad-i = TX(T), is 5 — d, the MF approximation 
breaks down when A^-i ~ m 5 ~ d . From this, and 
with the help of (0), one obtains the following 



criterion for the applicability of the MF scaling: 

In the special case of d = 4 one finds: m ^> 
T/[N\n(A/T)]. 

The Ginzburg criterion (|^) says that for masses 
m inside the window T ^> m ^> T/N x the 
MF scaling holds, while for smaller masses m <C 
T/N x (i.e., closer to the transition) the non- 
trivial d — 1 Ising scaling sets in. We see that 
the size of this latter, non-trivial critical region is 
suppressed at large N. 

3. Lattice 

The strategy we use to confirm our analytical 
results is the following. First, we show that the 
behavior at sufficiently large correlation length is 
given by the universality arguments, which in this 
case predict Ising Z(2) critical exponents. Sec- 
ond, we identify the boundary of the MF region, 
test the scaling and extract the exponent x. 

We discretize the d = 3 theory ([|) on a cubic 
L 2 S x L t lattice in the following standard way |J: 

N/2 

x 2—1 x,y 

x (x,x) 

3.1. The FSS results 

In order to study the critical behavior on lat- 
tices available to us we use the finite size scaling 
(FSS) method. Unlike previous bulk measure- 
ments of the critical exponents, when one has to 
keep the correlation length much smaller than the 
size of the box, in the FSS we let the correlation 
length saturate at the box size. This allows us to 
reach the Ginzburg region. 

The results of the FSS analysis for L t = 6 lat- 
tices are summarized in Table B. We see that all 
measured exponents show that the scaling very 
close to criticality is that of the Ising model in 2 
dimensions rather than the MF one. 

3.2. The crossover exponent x 

A straightforward way to find the value of x 
is to study the dependence of the order param- 
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Table 1 

Summary of FSS results and comparison with 
Ising model and MF scaling behavior. 



Exponents 


FSS 


Z(2) 


MF 


V 


1.00(3) 


1 


0.5 


0m /v 


0.12(6) 


0.125 


1 


7/1/ 


1.66(9) 


1.75 


2 



eter, £, on (3. Since £ vanishes at the critical 
point, it can be thought of as a measure of the 
distance from the criticality. We expect that for 
sufficiently small S, i.e., close to /3 C , this depen- 
dence should be given by a power-law scaling of 
the 2d Ising model: S ~ (const — /3) 1 / 8 . For 
larger £, further away from the criticality, the 
MF scaling holds: £ ~ (const' - /3) 1/2 . For even 
larger £ we should see the scaling corresponding 
to the fixed point of the 3d Gross-Neveu model 
§: £~ (const"-/?) 1 . In our simulations we can 
clearly resolve the MF region with the crossover 
towards the 2d-Ising region (Figure |l|) . 

In order to find the exponent x we use the re- 
lation: £ = 1/(21!) valid in the MF region. Using 
the histogram reweighting method and monitor- 
ing the quality of the linear fit as we add or sub- 
tract points we find for the boundary of the MF 
region: £ M f = 0.377(11), for N = 4, £ M f = 
0.213(4) for N = 12, and £ M f = 0.168(20) for 
N = 24. It is clear that the non-trivial 2d- 
Ising region is squeezed as N increases. The fit 
Smf = const ■ N~ x gives x — 0.51(3) which is in 
agreement with the analytical prediction x = 0.5 
ofeq. (§. 

4. Discussion and conclusions 

We have shown that a suppression of the crit- 
ical region occurs in theories with large number 
of fermion species. The corresponding crossover 
exponent x can be also determined (||). 

The role of the fermions is to screen the ef- 
fective self-coupling of the scalar field, A. The 
strength of this effect depends on two factors: (i) 
large N, and (ii) large window of scales between 
the cutoff of the effective theory, A, and the tem- 
perature, T. In QCD there is almost an order 
of magnitude window between the scale of the 
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Figure 1. Order parameter squared vs. f3 for 
lattice theories with N = 4, 12, 24 (L t = 6). The 
straight lines are the fits to the data in the MF 
regions. 



spontaneous symmetry breaking A ~ 1 GeV, and 
T c « 160 MeV, which is presumably sufficient to 
drive the effective self-coupling of the scalar field 
to its infrared fixed point value at the scale of T c . 
How small this value is then depends on the num- 
ber of the fermions (the condition (i)). The value 
N c — 3, though not very large, can be consid- 
ered large in some cases. One can expect, there- 
fore, that in QCD the width of the critical region, 
where actual 0(4) exponents are observed is sup- 
pressed. Further analysis is required to make a 
quantitative prediction. 
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